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We study the phase diagram of two models of spin- 1/2 antiferromagnets composed of corner- 
sharing tetrahedra, the basis of the pyrochlore structure. Primarily, we focus on the Heisenberg 
antiferromaget on the checkerboard lattice (also called the planar pyrochlore and crossed-chains 
model). This model has an anisotropic limit, when the dimensionless ratio of two exchange con- 
stants, J x /J <C 1, in which it consists of one-dimensional spin chains coupled weakly together in a 
frustrated fashion. Using recently developed techniques combining renormalization group ideas and 
one-dimensional bosonization and current algebra methods, we show that in this limit the model 
enters a crossed diraer state with two-fold spontaneous symmetry breaking but no magnetic or- 
der. We complement this result by an approximate "quadrumer triplet boson" calculation, which 
qualitatively captures the physics of the "plaquette valence bond solid" state believed to obtain 
for J x /J ~ 1. Using these known points in parameter space, the instabilities pointed to by the 
quadrumer boson calculation, and the simple limit J x /J > 1, we construct a few candidate global 
phase diagrams for the model, and discuss the nature of the quantum phase transitions contained 
therein. Finally, we apply our quasi-one-dimensional techniques to an anisotropic limit of the three- 
dimensional pyrochlore antiferromagnet, an approximate model for magnetism in GeCu2 04. A 
crossed dimer state is predicted here as well. 



T3 

C 

O 

i o i 

(N 
> 

On 
(N 

m 
O 
in 
o 

c3 



O 

o 



X 



I. INTRODUCTION 

We consider one of the most frustrated two- 
dimensional (2D) antiferromagnets: the checkerboard an- 
tiferromagnet, also known as the planar pyrochlore and 
the crossed-chains model (CCM). As the name suggests, 
this model is motivated by the three-dimensional (3D) 
pyrochlore materials. The 2D model is obtained by a 
projection of 3D corner-sharing lattice of tetrahedra on 
a 2D plane. This projection maps a four-spin tetrahe- 
dron onto a four-spin square with additional links (anti- 
ferromagnetic exchanges) along the diagonals. The struc- 
ture obtained in this way, depicted in Fig. ^ preserves 
the corner-sharing arrangement of crossed squares, typ- 
ical of the original 3D pyrochlore lattice, but destroys 
the symmetry between bonds of the tetrahedron: in two 
dimensions, the horizontal and vertical bonds are not 
equivalent to diagonal ones. This lowering of symme- 
try suggests consideration of extended 2D models with 
the checkerboard structure where exchange interactions 
on horizontal/ vertical and diagonal bonds take on dif- 
ferent values. Among these, the quasi-one-dimensional 
limit, in which exchange along horizontal and vertical 
directions J is much stronger than that along diagonal 
directions J x , is of special interest because it involves 
competition between strong quantum fluctuations, typ- 
ical for one-dimensional (ID) spin chains, and equally 
strong geometric frustration encoded in the structure of 
the crossed-chains lattice. 

The resulting checkerboard antiferromagnet has been 
analyzed by a variety of techniques along several com- 
plimentary "directions" in the parameter space: semi- 
classical analysis in the limit of large spin S 




FIG. 1: (Color online) Heisenberg antiferromagnet on the 
checkerboard lattice, viewed as coupled spin chains. Horizon- 
tal (vertical) spin chains run along the x (y) axis. Spins be- 
longing to the horizontal (vertical) chains are shown as green 
(red) filled circles. The intra-chain exchange (thick lines) is 
J, and the inter-chain exchange (diagonal thin lines) is J x • 



large- N expansion^i& easy-axis generalization 
(of the 3D modcl)£ and a quasi-lD (J x /J -C 1) 
approach^ In parallel with analytic approaches, the 
model was investigated numerically via exact diagonal- 
ization studies&ifiiAI and cluster-based strong-coupling 
expansion techniques i&l&l^ The present paper comple- 
ments these approaches by combining a controlled anal- 
ysis of the quasi-lD limit with general arguments to pin 
down limits of the phase diagram and postulate a likely 
global phase structure of the model. 
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We begin by expounding the more general context 
of the problem. One of the central theoretical motiva- 
tions behind the study of frustrated quantum magnets 
is the hope that, when magnetic ordering is suppressed 
by frustration, more novel types of order or even critical- 
ity may emerge. Phenomenological approaches suggest 
possible interesting quantum phases exhibiting "valence 
bond solid" (VBS) order, in which spins pair into singlets 
that are spontaneously localized on specific bonds, break- 
ing lattice symmetries. More exotically, such approaches 
suggest the possibility of phases with "topological or- 
der" , in which spins fluctuate quantum mechanically in 
a liquid-like state with however subtle topological prop- 
erties and often excitations with anomalous (e.g., frac- 
tional) quantum numbers. More recent predictions from 
such theories also include "deconfined" quantum critical 
points and phases in which several types of quasi-long- 
range (power-law) orders coexist unconnected by micro- 
scopic symmetries. 

Unfortunately, these types of phenomenological meth- 
ods do not give precise guidance as to the specific models 
in which such quantum orders appear, and attempts to 
find them in realistic microscopic Hamiltonians have met 
with at best limited success. The one specific context in 
which examples of all the above phenomena arc, however, 
known to occur is in one-dimensional spin chains. More- 
over, the theoretical and microscopic understanding of 
such spin models is vastly more complete than in two or 
three dimensions. A natural hunting ground for the ex- 
otic phenomenology described above would hence seem to 
lie in spin models consisting of chains weakly coupled into 
two or three dimensional arrays. A recently gained un- 
derstanding of the crucial role of nominally irrelevant op- 
erators and fluctuation-generated interactions in describ- 
ing frustrated quasi-lD magnetic systems^ described 
below, brings the hunt to (some degree of) fruition. 

In this paper, as in a previous work^ we follow this 
approach, taking as the weakly coupled units in ques- 
tion 5 = 1/2 Heisenberg nearest-neighbor antifcrromag- 
netic chains (other further-neighbor interactions along 
each chain may be included, provided they are not overly 
strong). A cause for hope is that such a ID chain is well 
known to exhibit a critical ground state with power-law 
correlations of various types. One prominent type of cor- 
relation in such a chain is antifcrromagnctic, specifically: 

(S(n)-S(n'))~ t^ +■•• , (1) 
|n — n | 

where n is the coordinate along the chain, and the brack- 
ets indicate a ground state expectation value. The omit- 
ted ■ • ■ terms decay much faster (~ lj\n — n'\ 2 or faster) 
than the dominant slowly-decaying antifcrromagnctic one 
shown here (we have also for simplicity neglected an 
unimportant multiplicative logarithmic correction to this 
term). The dominance of antifcrromagnctic correlations 
in the two-spin correlation function often leads to the 
misconception that a good picture of the ground state of 



the Id Heisenberg chain is that of fluctuating local an- 
tiferromagnetic order, i.e., a magnet in which spins are 
locally Neel ordered but the quantization axis fluctuates 
in space and time. Such a picture is in fact incomplete. 
This becomes clear upon considering the fluctuation of 
the local bond energy or dimerization, 

B(n) = S{n) ■ S{n + 1) - (S(n) ■ S(n + 1)). (2) 
One finds that its (staggered) correlations, 

Wn)B{n'))~t£^ + --., (3) 

have precisely the same slow power-law decay (again, up 
to a multiplicative logarithmic correction) as the anti- 
ferromagnetic ones in Eq. above! Further examina- 
tion of other correlators reveals no additional power-law 
correlations with competitive slow decay. Thus the ID 
Heisenberg chain should be thought of as consisting of 
locally fluctuating antiferromagnetic and valence bond 
solid order of comparable strength. 

With this understanding, it is natural to expect that 
weakly-coupled arrays of such chains might be pushed by 
the inter-chain coupling into magnetically ordered, dimcr 
ordered, or perhaps critical states, if this coupling favors 
the intrinsic antiferromagnetic or VBS ordering tendency, 
or fosters their balanced competition, respectively. While 
we believe this reasoning to be essentially correct, for 
many years, the richness of such possible behaviors went 
unrealized in the literature. This is because if the spin 
chains arc linked by magnetic two-spin Heisenberg inter- 
actions, these couple primarily to the antiferromagnetic 
fluctuations within the chains, and not to the VBS ones. 
Hence, for such a case, the problem of Heisenberg spin 
chains coupled by weak inter-chain interactions is rather 
well understood. With non-frustrated transverse (with 
respect to the chain direction) couplings, both renormal- 
ization groupie and self-consistent mean-field analysis^ 
predict an instability towards classical long-range ordered 
phase characterized by a non-zero expectation value of 
the spin (S r ) ^ 0. This instability follows from the cor- 
relations in Eq. QJ, which, loosely speaking, make the 
spin chain highly susceptible to magnetic ordering. 

More recently, it was recognized that the situation be- 
comes more interesting and less clear-cut when the inter- 
chain interaction is strongly frustrated, as is the case 
for the crossed-chains model we investigate here. The 
effect of frustration is to reduce and, ultimately, nul- 
lify the effective inter-chain magnetic field experienced 
by spins of the chain due to transverse inter-chain ex- 
change interactions (due to cancellations between contri- 
butions from spins whose local fluctuating orientations, 
according to Eq. are antiparallcl) . With no effec- 
tive external field present, the classical ordering insta- 
bility is naturally absent, resulting in (almost) decou- 
pled behavior of distinct spin chains. In formal calcula- 
tions embodying this physical picture, the weak residual 
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inter-chain interaction which does not cancel with pre- 
dominantly antiferromagnetic correlations appears to be 
described by the scalar product of conserved spin cur- 
rents from the chains involved. This observation led to 
the proposal that, as a result, the system of such cou- 
pled chains forms a liquid-like ground state with frac- 
tionalizcd spin excitations (spinons). The systems con- 
sidered included a frustrated spin ladder its 2D ex- 
tension, i.e., the spatially-anisotropic frustrated square 
lattice antiferromagnet^ and the crossed-chains models 

As shown in Ref . 0, in the former two cases these con- 
clusions arc in fact incorrect, due to the neglect of the 
VBS correlations, Eq. ©, equally as inherent as the an- 
tiferromagnetic ones to the Heisenberg chain. Although 
the microscopic magnetic exchange between spins on dif- 
ferent chains does not directly couple to VBS fluctua- 
tions, such a dimer coupling between (certain pairs of) 
chains is inevitably generated by the weak residual mag- 
netic interactions remaining after the dominant antifer- 
romagnetic cancellation. A careful analysis of the types 
of such dimer couplings allowed by symmetry and the 
detailed mechanism of their generation are crucial in de- 
termining the fate of the spin system and the strength of 
any ordering tendency. 

Technically, this analysis can be accomplished in a con- 
trolled fashion using powerful field-theoretical methods 
borrowed from ID physics. The point, made in Ref. Il5l 
is that no fine-tuning of the two-spin inter-chain ex- 
change interaction can make the low-energy field the- 
ory exactly of current-current type. Some higher-order 
derivative terms (typically involving spatial derivatives 
of the staggered magnetization field) are bound to be 
present (as getting rid of all of them to all orders would 
require tuning infinite number of inter-chain couplings to 
zero). Such derivative terms are commonly neglected on 
the grounds of their irrelevance with respect to the Lut- 
tingcr liquid fixed point of the independent spin chain. 
However, the quasi-lD problem is not the same as the 
purely ID one. Instead of disregarding irrelevant high- 
derivatives terms from the outset, one has to consider 
if they, in combination with the leading current-current 
term, can produce quantum corrections to the relevant 
inter-chain couplings. This indeed occurs both in the 
models of Ref. [jjj, and, as we will see, in the crossed- 
chains model studied here. 

In the present paper we extend the analysis of Ref. ^| 
to the CCM and show that previous claim of the sliding 
Luttinger liquid ground stated is not correct. Instead, 
similarly to the spatially-anisotropic square lattice model 
discussed above, the ground state is of spontaneously 
dimerized type, albeit with staggered ordering of dimers 
on parallel chains. The resulting configuration, shown in 
Fig. can be described as a crosscd-dimer one. 

The paper is organized as follows. Section ITT1 describes 
the Hamiltonian of the CCM model, its lattice symme- 
tries, and the passage to the field-theoretical description 
of the low-energy degrees of freedom and the operator 
product algebra they form. Section IlIII describes pertur- 




FIG. 2: (Color online) Crossed-dimer dimcrization pattern. 
"Strong" bonds (ones where e > 0) on horizontal (vertical) 
chains are shown in green (red). As before, spins on horizontal 
(vertical) chains are denoted by green (red) circles. 



bative analysis of the model in the one-dimensional limit 
of weakly coupled chains, J x /J <C 1. It contains key 
technical details of our work and explains the mechanism 
by which the crossed-dimer phase is stabilized. The limit 
of the fully two-dimensional model (J x ~ J), the pla- 
nar pyrochlore antifcrromagnct, is analyzed within the 
plaquette-operator mean-field approximation in Sec. II VI 
This is followed by Sec. [Vj which summarizes the pre- 
ceding material in terms of two possible scenarios for the 
global zero-temperature phase diagram of the checker- 
board antiferromagnet. There we present phenomeno- 
logical symmetry-based analyses of the quantum phase 
transitions between various phases of the model (and also 
point out an interesting connection with the recent de- 
confined quantum critical point idea) . Section IV1I de- 
scribes a three-dimensional extension of our model, the 
quasi-one-dimensional pyrochlore antiferromagnet, and 
its possible relevance to the experiments on GeCu204 
and ZnV204. Our main points are briefly summarized 
in Sec. IVIII Two Appendices contain important techni- 
cal details of the fermionic formulation of the low-energy 
sector of the 5=1/2 isotropic Heisenberg chain. 



II. FROM LATTICE TO CONTINUUM FIELD 
THEORY 

A. Lattice model and symmetries 

The Hamiltonian of the system H describes a collection 
of horizontal (Hh) and vertical (H v ) Heisenberg chains 
interacting with each other via the inter-chain interaction 
V: 

H = H + V = H h + H v + V. (4) 
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Spins (S = 1/2) are located at the sites of the checker- 
board (crossed-chains) lattice shown in Fig.Q] The cross- 
ings of the lattice have integer coordinates (n, m), so the 
sites of horizontal chains have half-integer cc-coordinatcs 
n + \ and integer y-coordinate m, while sites of the ver- 
tical chains are described by (n, to + i) pairs. With this 
convention the Hamiltonian of horizontal chains reads 

H h = jJ2S h (n- 1/2, to) • S h (n + 1/2, to). (5) 



Similarly, H v is given by 

H v = J^S v (n,m- 1/2) • S v (n,m + 1/2). (6) 

With local uniform magnetization defined by 

s h (n,m) = S h (n - 1/2, to) + S h (n + 1/2, m), (7a) 
s v (n,m) ~ S v (n, m — 1/2) + S v (n, to + 1/2), (7b) 



with the action (in 1 + 1 dimensions) 22 23 



the inter-chain interaction reads 

V = J x s h (n, to) • s„(n,i 



(8) 



and is characterized by the inter-chain exchange J x > 
which is much smaller than the in-chain antiferromag- 
netic exchange J > 0. We note that J x is the nearest- 
neighbor exchange on the checkerboard lattice while J is 
the next-nearest-neighbor exchange interaction. 

The space group symmetries of H, Eq. can now be 
summarized. The translational subgroup is generated by 
unit translation along the horizontal chains Th and that 
along the vertical chains T v . The remainder is generated 
by 7r/2 rotations about a crossing, and reflections about 
e.g., a vertical line through either a site or midpoint of 
a bond of a horizontal chain. We denote these two op- 
erations "site parity" P s h and "link parity" Plk-, respec- 
tively. As these are microscopic lattice symmetries, they 
will be preserved by any renormalization group trans- 
formation. Observe that Pl is a product of two other 
operations: P Lh = P sh ° T h . 



B. Continuum field theory and scaling operators 

The limit J x <^ J allows us to approach the problem 
from one-dimensional perspective: we treat V as a per- 
turbation and ask whether it can destabilize the critical 
ground state of the independent (decoupled) spin chains. 
The smallness of the J x /J ratio allows us to take the 
continuum limit along every chain involved. As men- 
tioned in the Introduction, a single Heisenberg chain is 
described in the continuum limit (i.e., at low energies) 
by a universal critical theory, with a variety of power-law 
correlations. Formally, this is most compactly described 
as the Wess-Zumino-Witten (WZW) SU(2) 1 theory 2 !, 



,5' 
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WZW 



d 2 x Tr d^g^d^g 
~lk f<f x ^^Jd^d„gg^d x g. (9) 

Here g is an SU(2) matrix. The coordinate xq = vt (v is 
the spin velocity and r is imaginary time) and x\ = x, the 
coordinate along the chain, and d 3 x is defined by extend- 
ing this 2D space into a three-dimensional hemisphere 
x-i < 0, the boundary of which is the (compactified) 
physical 2D plane (xo,xi), and analytically continuing 
g(xo,xi) into this hemisphere such that g{x$, x\, x-i — > 
— oo) — > 1 and g{xa, x\, x-i = 0) = g(xo,xi). This for- 
mal action is not of very much direct practical use, but 
serves to illustrate the underlying degrees of freedom of 
the critical theory. All operators in the WZW SU(2)i 
theory can be constructed from g. Corresponding to the 
two dominant power-law correlations in Eqs. and 
there are two scaling operators^ 2 * 2 ^ 



N ~ —iTrga, 
e - Trg. 



(10) 
(11) 



Here a is the vector of Pauli matrices, and the ~ indicates 
that the proportionality between these fields and the 
physical staggered magnetization/dimerization involves 
a cut-off dependent factor. The operator N represents 
the local staggered magnetization, while e represents the 
local staggered dimerization (it is the continuum version 
of the bond operator in J2J). There are also subdomi- 
nant power-law correlations arising from fluctuations of 
the chiral SU (2) currents, 



Jr 
Jl 



1 

4tt 
1 

4tt 



Ttg^dga, 
Tr gdg^a, 



(12) 
(13) 



with d = (d - idx)/2, and d = (d + idi)/2. Phys- 
ically, the operator J = Jr + Jl represents the local 
uniform magnetization, while u(Jr — Jl) represents the 
local magnetization (spin transport) current. 

All the low-energy power-law correlations of the weakly 
coupled Heisenberg chains can be exposed by decompos- 
ing lattice operators into a set of the above continuum 
operators (and generally their derivatives, see below) for 
each chain. This, for example, leads to the following 
decomposition of the spin at a site n — 1/2 along the 
horizontal chain number to: 

S h (n - 1/2, m) = a [f h , m (x) + (-1)" N h ,m(x)] ■ (14) 

Here x = na (a is the lattice spacing) and J (respectively, 
N) represents uniform (respectively, staggered) part of 
the spin density. Similarly, for the vertical spin chains 
we have 

S v (n, to - 1/2) = a \f Vtn (y) + (-l) m N v , n (y)] , (15) 
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where y — ma. Notice that the continuum limit is taken 
only for the coordinate along the chain; the perpendic- 
ular one becomes an index m (respectively, n for ver- 
tical chains). The uniform spin magnetization J is the 
sum of the right ( Jr) and left ( Jl) moving components, 
J = Jr + Jl, and represents the conserved spin density 
(it is often referred to in the literature as the spin "cur- 
rent" , the term originating from the relativistic concept 
of space-time current, whose time component is the con- 
served density). Note that the staggered dimerization e 
does not appear in Eq. i|14|) ; in fact, it cannot appear in 
the decomposition of any single spin operator since it is 
not a vector under SU(2). As discussed in the Introduc- 
tion, for this reason dimer order does not appear likely 
in weakly coupled Heisenberg chains with unfrustrated 
inter-chain couplings. 

The action of the microscopic space group symmetries 
(described above) upon the continuum scaling operators 
will be crucial in the following. These are rather clear on 
physical grounds^: 

Translation: 



J, N 



-N, 



T : J 
Site parity: 

P S :J^J, N^N, e->-e. 
Link parity: Using Pl = P s ° T we find 

P L : J -> J, N -> -N, e -> e. 



(16) 
(17) 
(18) 



We will see at the end of this section that this symmetry 
is responsible for the absence of N v ■ Nh terms in the 
Hamiltonian of the problem. 

Because of somewhat non-intuitive point-splitting 
identities, the WZW model can be written in Hamilto- 
nian form (known as the Sugawara form) in terms of the 
spin currents. For a single chain, one has 



2ttv 

Hwzw = 



dx 



Jr{x) ■ Jr{x) + J L (x) ■ Jl(x) 



(19) 

Applied to the set of horizontal chains (labelled by m), 
the lattice Hamiltonian Hh, Eq. JjJJ, transforms into 



dx Jh,m,R{x) ■ Jh 

,iti,r{x) 

+ Jh,m.,L{x) ■ Jh,m.,L{x) 
+ 9bsJh, m ,R( X ) ■ Jh,rn,L{x) ■ (20) 



Here v = \Ja is the spin velocity. Note again that 
J£ m R i L (x) depends on position x — na along the chain 
direction whereas its y = ma coordinate dependence only 
shows up via the (horizontal) chain index m. We have ac- 
tually included in Eq. (|2U[) a correction (proportional to 
gbs) to the WZW model, which is present in the Heisen- 
berg chain but is marginally irrelevant in the situation 



under consideration. For this reason, it may be safely 
neglected in what follows. Similarly 



2tyv C r -> -> 

~3~X! / dy J v,n,R(y) ■ Jv,n,R(y) 
n 

+ Jv,n,L(y) ■ Jv,n,h{y) 

+ gbsJv,n,R.{y) ' Jv,n,L(y) ■ (21) 

C. Decomposition of the full lattice model 



Now we are ready to express the inter-chain perturba- 
tion Eq. JHJl in terms of low-energy modes J and N. We 
begin by analyzing the sum of two neighboring spins on 
the same (say, horizontal) chain, 

s h {n,m) = §h(n— 1/2, m) +Sh(n+ 1/2, m) 

= a [2J h , m (x) - {-l) n ad x N h , m (x)\ . (22) 

For the reasons to be explained in detail below, we 
have retained the next-to-leading irrelevant contribution 
(d x N) in this expression. Similar decomposition is done 
for the sum of two spins on the crossing vertical chain. 
The interchain interaction V thus reads 



V 



^2 {9jjJh,m( x ) ' Jv,n(y) 



- 9nj (-l) n d x N h , m {x) ■ Jv,n{y) 

+ (-l) m Jh, m (x)-d y N Vtn {y) 

+ g n n(-l) n+m d x N h . m (x) ■ d y N Vtn (y)} , (23) 

where, as before, x — na, y — ma and the following cou- 
plings arc introduced to shorten notations: 



f-Jjj — 4Jx a , g,ij — 2J x a , g nn — J^a . 



(24) 



It is important to observe that Eq. (|23|l does not con- 
tain Nh ■ N v type of terms, which are forbidden by the 
symmetry of the checkerboard lattice. For example, re- 
flection with respect to the vertical chain changes sign of 
Nh, {Pl ■ Nh — » —Nh), while leaving N v invariant, see 
Eq. (|18fl . This reflects strong frustration of the model 
under study, as discussed in the Introduction. Observe 
also that any pair of horizontal and vertical chains cross 
only once, which makes Eq. (12.'.>[| local in space. This re- 
quires us to think carefully about the short-distance reg- 
ularization of the low-energy theory defined by Eq. 12011 , 
Eq. (|21|) . and Eq. Il2.'j[| — the corresponding analysis is 
described in the next Section. 



D. Operator product expansion 

Various perturbations to the WZW model Eq. i|19|) 
[such as the intra-chain backscattering gi, s in Eq. I|20|) 
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and Eq. (|21() . and the inter-chain V, Eq. (|23|l ] are most 
conveniently analyzed with the help of operator product 
expansions (OPE). These are operator identities that are 
derived by applying Wick's theorem to a correlation func- 
tion of a pair of operators at nearby points, say, (x, r) and 
(0, 0) — several of the examples below are worked out in 
Appendix [B| sec also Appendix A of Ref. |2f| for more 
examples. The OPE below are valid for operators from 
the same chain, and, to lighten expressions, we suppress 
chain indices here. 

The spin currents Jr/l obey the following chiral OPEs, 
which arc frequently used in the literature^ [these, for 
example, are used to derive the renormalization-group 
flow of gbs term in Eqs. f2T)|t and QI) ]: 



J£Or,r)4(0) 



r<i b 



and 



J a L (x,r)J b L (0) 



8tt 2 v 2 (t — ix/v + aa T ) 2 

^K/|(o) 

2t:v(t — ix/v + aa T ) 



gab 

8it 2 v 2 (t + ix/v + aa T ) 2 
ie abc Jj{0) 
2itv(t + ix/v + aa T ) ' 



(25) 



(26) 



where, as explained in Appendix IA 21 a = a/v is the 
short-time cutoff of the theory and a T = sign(r). 

Being a conserved current, J a is also a generator of 
rotations. Thus the OPE of J a and JV" should be non- 
trivial. In fact, this one is the most important OPE for 
the subsequent analysis (see Appendix iBl for the deriva- 
tion), 



J%(x,t)N"(0) 
J£(x,t)A"(0) 



ie abc N c {0) - iS ab e(0) 



47rw(r — ix/v + aa T ) 
ie abc N c (0) + iS ab e(0) 
At:v(t + ix/v + aa T ) 



(27) 



Finally, fusing spin current with dimcrization e gives 
back the staggered magnetization 



J£Or,r)e(0) = 
J2(x,r)e(0) = 



iN a {0) 



■At:v(t — ix/v + aa T ) ' 

-iN a (0) 
4:Ttv(t + ix/v + aa T ) 



(28) 



Observe that Eqs. I|25|) - (|28J) form a closed operator 
algebra — this is the key technical reason behind the gen- 
eration of the inter-chain interaction of staggered mag- 
netizations in frustrated spin chains models, see Ref. [l5| 
and Sec. UTTBl below. 



III. LOW-ENERGY HAMILTONIAN 

The spatially-anisotropic J1-J2 modeU^ has taught us 
that keeping track of the nominally irrelevant terms is 



crucial for a correct solution of the problem. In this sec- 
tion, we extend this line of thinking to the crossed-chains 
model and demonstrate that indeed irrelevant terms pro- 
duce symmetry-allowed relevant ones in a simple pertur- 
bation theory. 



A. Symmetry analysis 

Before proceeding with microscopic calculations, it is 
instructive to write down most general form of the inter- 
chain Hamiltonian SV which is allowed by symmetries of 
the crossed-chains lattice. The reason to do so is that, 
while many such terms will be absent in a naive contin- 
uum limit of the original spin model, those which are "ac- 
cidentally" missing (i.e. not prohibited by any symme- 
try) may be expected to be generated as a "quantum cor- 
rection" (i.e. through a RG transformation) when naively 
irrelevant terms are taken into account. The necessary 
complete set of space group generators for this analysis, 
Th, P s h, Plh, and was introduced in Sec. Ill Al 

Naturally (as in any field theory), there are an infi- 
nite number of possible interactions, and since there are 
additionally an infinite number of chains, the multitude 
of potential terms is compounded. Physically, however, 
"pairwise" interactions involving fields on only two chains 
at a time are expected to be most important (interactions 
involving more chains simultaneously can be shown to oc- 
cur only in higher order in J x /J). Such an inter-chain 
Hamiltonian naturally splits into the sum of 5V X , which 
describes interactions between two crossing chains, and 
<5Vj|, which includes interactions between parallel chains, 
SV = SV X + SV\\ . Within these chain-pair interactions, 
we narrow the search by considering the "most relevant" 
possibilities (ones involving the smallest number of the 
smallest scaling dimension primary fields N and e and 
no derivatives). Since we are perturbing the decoupled- 
chain system, the appropriate sense of "relevant" is that 
of the decoupled 1+1-dimensional critical theories. We 
find 



SV* 



a x {-l) n + m e Km {na) e v , n (ma) (29) 



and 



SVn 



^ \a 2 {l)N v , m {na) ■ Nu,m+i(na) 



n.m.l v—h.v 



+ a 3 (l) (na) e v , m +i(na) ■ (30) 

We note that in Eq. (|30|l . an interaction is possible be- 
tween parallel chains an arbitrary distance I apart. From 
the point of view of the decoupled-chain fixed point, there 
is no notion (or effect in RG rescaling) of "distance" be- 
tween chains, so all such terms are equally "relevant" in 
this point of view. One would expect, however, these 
terms (i.e. 0,2(1), ^3(0) to decay in magnitude with in- 
creasing I. 
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It is straightforward to check that these terms and 
only these terms satisfy the symmetry requirements of 
the checkerboard lattice. First, the invariance of 6V\\ is 
easy to establish, as it involves pairs of operators e and 
N from like chains (i.e., horizontal- horizontal or vertical- 
vertical). These transform identically under all opera- 
tions, and invariance is trivially shown. 

The crossed-chains term, 8V X , is more involved. We 
sketch the arguments for its invariance. Rotation by ir/2 
about a crossing is manifest, as the fields in Eq. 1|29|) 
are drawn from a single such crossing. Unit translation 
along the x-direction makes eh —* —£h while e v is obvi- 
ously not affected. However, (— l)"+ m also changes its 
sign, (-1)"+™ -> go that T h (6V x ) = SV X . 

Reflection with respect to a site on a horizontal chain 
P s h preserves e v but does change sign of dimcrization on 
every horizontal chain: P s h(^h) = — e h- But at the same 
time P s h interchanges even and odd vertical chains, i.e., 
Psh((-iy i+m ) = -(-l)" +m . Thus P sh (8V x ) = SV X . 
Link parity Plh is simple since every e is even under it. 
Moreover, since Plyi is nothing but reflection with re- 
spect to, say, the vertical chain number n, the vertical 
chain with index n + 1 then transforms into that with 
index n — 1, etc. Hence, even and odd vertical chains 
are not interchanged by Plh, and (— \) n + rn — > (— l)"+ m , 
showing the invariance under this final generator. Notice 
that the staggering factor (— l) n + m plays a very impor- 
tant role in this consideration - its presence makes the 
local interaction of staggered dimcrizations possible. 

One could wonder if SV X could similarly include a 
staggered product of magnetizations, (— l) n+m Nh ■ Nh, 
but this is prohibited by the Phh symmetry. We note 
that microscopically, such a term cannot be generated 
(see the following subsection for the mechanics of gen- 
eration of the allowed terms) as a consequence of the 
identity J a (x, r)e(x, t') = which follows from the OPE 
Eq. I|28|) . The only symmetry-allowed combination of 
TV's that can show up in SV X is (Nh • N v ) 2 . Such a term 
does arise in the large-S" "order-from-disorder" calcula- 
tions, see Ref. 0, but in the S = 1/2 microscopic model 
it has scaling dimension 2 and is thus deemed irrelevant. 
Moreover, one can derive, using abelian bosonization, 
the OPE of two TV fields at the same spatial point x: 
N a (x,r)N b (x,0) ~ ie afcc sign(V)J c (a;,0). This allows one 
to identify2& this biquadratic term with the dimension 
2 scalar product of two spin currents on crossing chains 
[that is, the gjj term in Eq. (J23J], (Nh • N v ) 2 — > J h ■ J v . 

Observe now that none of the symmetry-respecting 
terms in 5V X and 6V\i are present in the naive continuum 
limit of the theory Eq. I|23[) . Below we show that second- 
order perturbation theory in the inter-chain exchange 
J x generates SV X with coupling constant oi ~ J x /J- 
Similar arguments show that 02.3 ~ J x /J 3 . This is 
because in J 2 order one generates terms involving a 
product of derivatives of N fields on parallel chains, 
~ J x d x Nh ■ d x Nh- Once these are present, one can fol- 
low calculations in Ref. to find that both ai and 03 



terms in <5Vy are generated, but this happens only in the 
next, (J 2 ) 2 = J X ; order of the perturbation expansion. 
Since, as we show below in Sec. IIII CI the 5V X contri- 
bution is relevant, it is sufficient to keep only the lead- 
ing a\ type terms sublcading 02,3 ones arc too small 
(02,3/01 ~ J x /J 2 *C 1) to change the outcome. 



B. On the importance of irrelevant terms 

Here we describe microscopic calculation of SV X . We 
begin by expanding the action in powers of V Eq. J22J). 
This generates a number of terms of which the most im- 
portant ones involve products of spin currents and stag- 
gered magnetizations from the same chain (and with the 
same spatial coordinate — that is, all fields belong to the 
same crossing). Thus we pick out gjjg n n and cross-terms 
from g 2 -. These contributions can be written in the form 



EE(-r)" +m / dTdr'v x (x,n;y,m;T,T') 



a.b n,m 



where 



x—na^y—ma 

(31) 



v x =g* j [d x N^ m (x,T)Jl m (xy)J- n (y,r)d y N^ n (y,T') 

+ Jl m (x,T)d x Nl m (x,T')d y Nl n (y,T)jl n ( y y)\ 

+ 9n9nnJh, m (x, T)d x N h h m ( Xl t') 

xJl n (y,r)d y Nl n ( V y). (32) 

Now apply the OPE Eq. (|2*7Jl to the product of fields 
from the same chain and at the same spatial point x. For 
example, 

J^(x,T)d x N b (x,T') = lim d x ,J a (x,T)N b (x\T') 

x'—>x 

-i[ie abc N c (x, t') - iS ab e(x, t')] 



AtTV 2 (t — T 1 + UG r - T i 



(33) 



Similarly 

J a L (x,r)d x N b (x,r') 



i[ie abc N c (x, t') + iS ab e(x, t')] 
Attv 2 (t — t' + acr T _ r ') 5 



J34) 

We observe that the OPE of the full spin current J and 
N at the same spatial point does not contain staggered 
magnetization 

J a (x,r)d x N b (x,r') = [J a R (x, r) + J a L (x, r)]d x N b (x, r') 
= 5 ab e(x 1 T l ) 

2ttv 2 (t - t' + aa T - T ') 2 

This is a very important result, with which Eq. (|3 1|) can 
be brought into the surprisingly compact form: 

^2 - ^<5 Qb ^(-l)' l+m (2 5r 2 y+ffjj5rm ) 

a, b n 7 m 

, €h, m (na,T')e Vin (ma,T') 

x / alrdT - — — — —. (36) 

[27tw^(t — t + a<7 T - T iy\ z 
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The integral involved is obviously convergent 
1 2 2v 3 

dt 



{t + aat) 4 3a 3 3a 3 ' 



(37) 



Using Eq. I|24[l for the <?'s involved, we finally obtain 
the following fluctuation-generated correction to the low- 
energy effective action 



5S 



3 J? a 3 



J a f 

~5 — / dTy^(-l) n+m e h , m (na,T)e v>n (ma, 

-K A V J * — ' 



Denoting 



ZJla 2 



(38) 



(39) 



we have the following addition to the inter-chain Hamil- 
tonian V, Eq. (|23(l . to analyze [this is because Z = J e~ s 
and SS = J drSV] 

5V X =- Se «E(^) n+me "^(H£»,r,(H- (40) 

n,m 

Staggered dimerization e has scaling dimension 1/2, 
which means that it is as important for the chain physics 
as N is. In fact, up to logarithmic corrections, correla- 
tion functions of staggered dimerization and magnetiza- 
tion decay with the same power law, a; -1 , sec (I1I2II . 

This is also clear from the OPE Eqs. (H7J) and 
which show that N and e transform into each other un- 
der chiral rotations generated by Jr/l- Since any pair 
of horizontal and vertical chains has only one crossing, 
Eq. H4Uf) is a sum of local terms, each of which is marginal 
(space-time dimension = 1, and dimension of the product 
ehtv is 1 as well). However, we shall see below that this 
marginality is superficial an infinite number of marginal 
crossings add up to a relevant perturbation. 

C. Mean- field analysis of the effective inter-chain 
interaction Eq. H4UD 

From now on it is safe to omit derivative terms present 
in V, Eq. i|23[) : their role was to generate, as described 
in Sec. lIII 51 more relevant symmetry-allowed inter-chain 
interactions. With this in mind, we write down the renor- 
malized version of Eq. (|23|l . 

V = ^ \gjjJh, m (na) ■ J v ,n(ma) 

n.m 

- (-l) n+m g e ae htm (na)e V!n (ma) . (41) 

As discussed above, the first term originating from the 
naive continuum limit of Eq. © has scaling dimension 2 



while the second term, which is SV X generated by high- 
energy fluctuations, has scaling dimension 1. Thus we are 
allowed to discard the irrelevant current-current piece of 
V, Eq. H41fl . As a result, all that remains of the interchain 
interaction is given by SV X [Eq. (|40|l ] , V — ► 5 V x , which 
was not present in the naive continuum limit Eq. 123J) at 
all! We tackle it, in analogy with analysis of Ref. [Ij, by 
the chain mean-field approximation. The staggering fac- 
tor (—1)™+"' suggests a staggered dimer order on parallel 
chains. That is, we assume the pattern 

*h,M = e»,n(l/) = (-!)"<£>, (42) 

where (e) is a mean-field expectation value. The inter- 
chain coupling is then decoupled into a sum of indepen- 
dent single-chain Hamiltonians, 

SV X = -^(-l) m 5£ a( e )^e h ., m (.T = na) 

rn n 

- ^(-l)"o e a<e) Y, <*,»(» = ™«)- (43) 

n m 

Look on one of them, say, that of the horizontal chain 
with index m (which is fixed now). Now we can take 
continuum limit (J~) n f(na) — * a -1 J dxf(x)) and 

6V x (m) = -(-ir. 9e (e) J e(x) dx, (44) 

which can be easily analyzed along the lines of Ref. 0. 
Using abelian bosonization expression for the staggered 
dimerization 



e(x) = — cos[v / 27np(a;)], 



(45) 



where if is the spin boson field and A is a nonuniversal 
constant of order l^Si we arrive at the effective single- 
chain sine-Gordon action for the mth chain 

S(m) = J d 2 r(^(\7 v ) 2 -G cos V^tfY (46) 

The action S(m) is written in terms of dimcnsionlcss 
coordinates r = [x/a, vr/a) and the effective coupling 
constant G = A 2 g e (cos ^/2tt \p) / \tt 2 v) . The self-consistent 
equation for (cos y/2%(p) follows from the exact solution 2 ^ 
for the free energy F m of the sine-Gordon model Z m = 
J Dipexp[-S(m)] = exp(-F m ); 



(cos v 27T<p) 



dF d\nZ 



dG dG 3^3 
where the constant cq reads 

2r(l/6) /7iT(3/4)\ 2/3 



^G 1 / 3 , (47) 



Co 



/?rr(2/3) \2r(l/4) J 
Simple algebra gives 

^ /A 2 .g 



(cos w2-Kip) 



0.265316 



Jx 

J 



(48) 



(49) 
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where we have set A = 1. Hence the expectation value of 
the staggered dimerization is proportional to J x , 

(e) = — (cos%/2W> = 0.08444 1 . (50) 
na Ja 

The spin gap A is given by the mass to of the lightest 
breather in the sine-Gordon theory^ 



to = 2Msin(7r/6) = M, 



where M and G are related by 

r(i/4) 



G 
2~ 



7rr(3/4) 



A/ V^T(2/3)\ 3/2 
2r(l/6) J 



(51) 



(52) 



Thus M = cqG 2 / 3 , and finally the spin gap A is found 



as 



= M = 



4A 2 c, 3 



f ( —\ 



| -y 



(53) 

The resulting dimerization pattern is shown in Fig. 
An equivalent configuration is obtained by a global shift 
of crosses by one lattice spacing along either the x or 
y direction. It is worth pointing out that exactly such 
inter-dimer correlations - crossed-dimer ones - have been 
observed in the exact diagonalization study of finite CCM 
clusters, see Table II and Fig. 5 in Ref. 



idy ( 
11. 



IV. THE PLANAR PYROCHLORE: 
PLAQUETTE VBS AND ITS INSTABILITIES 

In the preceding sections we focused on the quasi-lD 
limit J x <C J, and established the existence of the spon- 
taneous long-range order of the crossed dimer configura- 
tion (Fig.EJ. In this section we will explore a different re- 
gion in the parameter space, where the nearest-neighbor 
coupling J x and next-nearest-neighbor exchange cou- 
pling J are nearly equal. Earlier numerical studies using 
exact diagonalizationifi and strong-coupling expansion 
tcchniquesi2ii&ii showed that the ground state at J = 
J x is a valence bond crystal with long-range quadrumcr 
order, shown in Fig. [3] Here we review a simple theoret- 
ical account of this plaquette VBS (P-VBS) state using 
the quadrumer-boson approximation^Si^Si^i and exam- 
ine its instabilities to other orders. This analysis, to- 
gether with the results in the preceding sections, will 
serve as a basis for our discussion on the global phase 
diagram of the CCM in the following section. Our sim- 
ple approach presented here is meant to give a qualita- 
tive picture; more quantitatively reliable numerical re- 
sults can be obtained, for example, by series expansion, 
as developed in Refs. and HI. 

In the following analysis it is more convenient to use a 
new coordinate system labelled by (J, k), rotated by 7r/4 
from the x and y axes; see Fig. |3J The quadrumcrized 
valence bond crystal breaks lattice translation symmetry, 




FIG. 3: (Color online) Quadrumerized checkerboard lattice 
with coordinates (j, k) shown. The plaquettes with blue cir- 
cles are quadrumerized. Each unit cell contains 4 spins. 



and each quadrumerized plaquette centered at (j, k) has 
four spins, Si (I = 1,2,3,4). From the outset we assume 
the breaking of translation symmetry and begin with the 
Hamiltonian for a single quadrumerized plaquette, 



Hp — J x ( Si ■ 02 + $2 • S3 + S3 ■ S4 + O4 • S\ 



Jx 
2 



Si + 02 + 5*3 + 5*4 
~ I So + S4 



Si + Si 



(54) 



The lowest-energy state of H p is a spin singlet with en- 
ergy — 2J X , which can be written as 



S T|0) 



1 



2\/3 



(lttu) + iutt> + ituT> + iim) 

- 21 T1T1> — 2| 1T1T », (55) 



where 1 01 02 03 04) denotes the state with Sf = ai- The 
first excited states are a triplet with energy — J x , 



tl\0) = i(| tm > + I TITT > - I TUT > - I ITTT )), (56a) 



4|o) = - 7 =(iuu)-imT)), 



(56b) 



* f -|o} = gGum + UTU}- IUTD- 1 Till})- (56c) 

The operators , t\, and t\ can be thought of as creation 
operators of hard-core bosons. 

As mentioned above, the ground state of the CCM is 
known to be a gapped P-VBS state at the planar py- 
rochlore point J = J x . As long as J w J x , we may 
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thus expect that a good approximation to the ground 
state should be obtained by direct product of the sin- 
glet states, Eq. weakly hybridized with the triplets, 
Eqs. (|56|) . Motivated by this observation, we employ the 
quadrumer-boson approximation22i22i2i in which we keep 
only the low-lying four states, singlet and triplet, in each 
quadrumcrizcd plaquctte, and discard the other higher- 
energy states. Now the boson operators are subject to 
the constraint 



sU + t\t + + t\t z + tit_ = 1. 



The plaquette Hamiltonian is then written as 



-2J X + J x {tU++tlt z +tU_ 



(57) 



(58) 



The spins Si can also be written in terms of the hard-core 
boson operators. The representations are found from ma- 
trix elements of the spin operators with the four states. 
After some algebra we find 



Sf 



s 



\ - *'-*_) + ^ (4* + *%) , (59a) 

L (tlt z - t\t_) tjS (f + s + sh_) , (59b) 

(tls + sh + ^) , (59c) 



1 

71 



tU. 



V6 

(zi) 

(-1)' 
V3 



where I = 1, 2, 3, 4. Assuming that the density of triplets 
is low in the P-VBS state, we keep only terms linear in 



tfj, and set s = = 1; 



Sf 



(-1)' 
V6 



(4 + *«) 



(60) 



where a — x,y, z and we have introduced 
1 , 1 



=(*+ + *-). *» 



(t+-t_). (61) 



With the coordinate system (j, k) in Fig. [31 the total 
Hamiltonian of the checkerboard antiferromagnet reads 

j,k 

+ Jx (»5j',M ' Sj-1^,4 + ^i,fc,2 • Sj t k-l,l 

3,k 

+ 5j,fe,3 ' "5j+l,fc,2 + Sj,k,4 ' ^7,fe-[-l,3 
■ Sj,k+l,3 + Sj,k,4 ' %fc+l,2 

+ Sj t k,3 ■ + Sj,k,4 ' "Sj+l,fc,2 



(62) 



where Sj,^; is the rth spin 5; in the quadrumerized 
plaquette centered at (j,k). With the approximation 



Eq. H60|) the Hamiltonian becomes 



p \ a=x,y,z / 



(63) 



where we have introduced the triplet boson field, 

Up) 
M-p)j 



(64) 



with the momentum p = (p\,P2) and the Fourier trans- 
form 



-i(jpi+kp 2 ) 



ta(j, k), 



(65) 



where Af is the number of quadrumerized plaquettes. 
The Hamiltonian matrix is given by 



f Jy.+e{p) e(p) 



where 



£ (.P) = o - ^x)(cOSpi + COSp 2 )- 



With the Bogoliubov transformation, 



ta(P) \ _ (coshOp sinh6» p -\ / b a (p) 
tU-P)J Unh% cosh9 p) W-P) 



where 



exp(-4dp) = 1 



Jx ' 



(66) 



(67) 



(68) 



(69) 



the Hamiltonian [Eq. (|63|) ] is diagonalized, 
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-E(p) + E($)J2bt(p)b a (P) 



(70) 

The energy dispersion of the triplet eigen mode b a (p) is 
given by 



E(p) 



Jx 



Jx + g(J - Jx)(cOSpi + COSp 2 ) 



1/2 



(71) 

The simple quadrumer-boson approximation described 
above reproduces the basic feature of the previous nu- 
merical studiesi 10 i 12 i 13 i 14 at J w J x the P-VBS state is a 
stable singlet ground state with a gap to excited states. 

The softening of the triplet mode Eq. lf7T|) tells us po- 
tential instabilities that the P-VBS state may have (see 
also Ref. B. First, from E(0) = [J x (8 J - 5J X )/3] 1/2 
we see that it becomes unstable at J x — > 8 J/5, when 
the bosons condense at p = 0. The resulting state 
has the long-range Nccl order, like the one realized in 
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FIG. 4: (Color online) The magnetically ordered state due to 
condensation of bosons b a (n,n). The solid circles represent, 
say, up spins, and the open circles represent down spins. The 
quadrumerized plaquettes in the nearby P-VBS phase are in- 
dicated by blue circles. 



the square-lattice limit J x 3> J. This can be eas- 
ily seen from Eq. 1|60[) . in which we may suppose that 
t\ + t a is a nonvanishing c-number upon condensation 
of bosons. In the present approximation, the transition 
from the plaquette singlet to the Neel-ordered state oc- 
curs at (J/Jx)ci = 5/8, which should be compared with 
the estimate (J/J x ) c i = 0.79 ~ 0.81 from a sophisti- 
cated strong-coupling expansion^ Second, the P-VBS 
state becomes unstable as J/J x — > (J/Jx)c2 = 11/8, 
at which bosons with momentum p = (jr, 7r) condense. 
Upon bose condensation the spins will have a magnetic 
long-range order, shown in Fig. with the spin config- 
uration TTIITTI4 along the diagonal directions and the 
Neel order along the horizontal and vertical chain direc- 
tions [as can be seen by replacing t< a + t a —* (— l) J+fc in 
Eq. 1)6(1 In fact, this is one of the magnetically ordered 
states which are shown to be stable at J 3> J x in the 1/5 
expansion; see Fig. 2(b) and Fig. 6 in Ref. 0. The soften- 
ing of the triplet mode at p = (jr, ]LL^ S a ^ so found in the 
strong-coupling expansion in Rcf. |lj with the estimated 
critical coupling (J/J x )c2 = 1-06 ~ 1.13. 

Finally, we comment on two defects in our approach. 
One is that the lattice translation symmetry is explicitly 
broken from the outset and cannot be restored within 
the theory. This means that the Neel ordered state at 
J x > 5 J/8 inevitably has the P-VBS order as well - 
i.e., it is a coexistence phase. This should probably be 
viewed as an artifact of the approach - see Sec.[V]for dis- 
cussion of this portion of the phase diagram. The other 
defect is that we have ignored interactions among the 
triplet bosons (as well as projected out the other higher- 




FIG. 5: (Color online) Global phase diagram of the CCM 
model in scenario I of Sec. IV Al Thick vertical lines with 
question marks indicate that the corresponding transition is, 
according to Landau theory reasoning, either first-order or 
occurs via an intermediate co-existence phase. 



energy states - singlet, triplet, and quintuplet - on each 
quadrumerized plaquette). In the crudest approxima- 
tion we adopted, the first excited states are a triplet ex- 
citation 5=1 with no energy dispersion at J x = J. 
The numerical studie3i2*iSiiii4 showed, however, that at 
J = J x the lowest-excited states are in the spin sin- 
glet sector (see, for example, Sec. V in Rcf. llOjh The 
cluster-based calculationsiSii^ indicate that these singlet 
excitations are bound states of two 5=1 excitations. 
To describe correctly the singlet excitations in terms of 
the quadrumer bosons, one needs to go beyond the linear 
approximation Eq. I|6(J|) and include interactions among 
the bosons. We do not try to do this here, but only point 
out that the dispersionless triplet mode is naturally sus- 
ceptible to forming a bound state. Away from the planar 
pyrochlore point J x = J, not much is known about the 
low-energy excitations in the 5 = sector; it is not well 
understood how the singlet energy levels in the spin gap 
change function of J x / J^a 



V. GLOBAL PHASE DIAGRAM OF THE CCM 

In this section we discuss the global zero-temperature 
phase diagram of the CCM as the control parameter 
J x /J is increased from to oo. Our analysis re- 
lies on three well-established facts. First, the ground 
state at J x ^ J obviously has the Neel order and is 
smoothly connected to the Neel ordered state of the an- 
tifcrromagnctic Heiscnbcrg model on the square lattice. 
Second, by now there is convincing numericaUS*i2iiii4 
and analytical^^ evidence for the P-VBS state at and 
around J = J x . Finally, as shown in Sec. II I II the 
ground state is spontaneously dimerized in the quasi- 
one-dimensional J x *C J limit as well, where long-range 
crossed-dimcr order (Fig. |3J) sets in. 

The first question we address here is how the two 
dimerized phases, plaquette VBS (P-VBS) and crossed- 
dimer VBS (CD- VBS), are connected. We propose the 
two complementary scenarios in subsections lV Al and IV Bl 
below. The nature of the transition between the P-VBS 
and the Neel states is discussed in Sec. IV CI 
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A. Scenario I: direct transition between the 
crossed-dimer and plaquette VBS 

The 'minimal' assumption is that the two quantum- 
disordered valence bond phases connect at some critical 
value Jx/J < 1- The validity of this assumption can only 
be verified by the exact calculation of the ground state 
energies of two dimerized phases in the whole interval 
< J x /J < 1 of interest. Such a calculation is obviously 
beyond our analytic approximations suited for J x / J — > 
(CD- VBS, Sec. E3 and J x /J -> 1 (P-VBS, Sec. llVl) 
limits. Instead, we take a phenomenological point of view 
here and assume that the ground state energies are such 
that a direct transition between the CD- VBS and P-VBS 
phases is possible. At least partial support to this point 
of view is provided by the exact diagonalization study 
of Sindzingre et ali^ which seems to indicate a single 
change in the ground state around J x /J ~ 0.8. 

The question then is if this transition between these 
two phases can be continuous. In general, this question 
is difficult to answer. Most formally, the rcnormaliza- 
tion group theory of continuous critical phenomena sets 
only some rather weak constraints on the existence of a 
continuous phase transition between any two phases "A" 
and "B" . In particular, it requires the existence of an 
abstract scale-invariant fixed point (critical field theory) 
with a single "relevant" symmetry-allowed operator in its 
spectrum, such that a positive/negative coefficient of this 
operator in the action takes the system into phase A/B. 
The critical fixed point theory itself must clearly have 
higher symmetry than either phase A or B, but no a pri- 
ori restriction is placed on the relation of the symmetries 
of phase A to those of phase B. 

A conventional - and more stringent - "criterion" for 
the existence of a continuous transition is based on the 
specific realization of the critical theory provided by a 
Landau-Ginzburg- Wilson (LGW) action written in terms 
of order parameters. More physically, LGW theory per- 
mits a continuous transition by the "condensation" of 
some "soft mode" of phase A, which transforms non- 
trivially under the symmetry group of A. The conden- 
sation of this soft-mode order parameter then leads to a 
lowering of symmetry (since by assumption the condensa- 
tion breaks some symmetries that it transforms under) in 
phase B. A necessary criterion for a LGW-allowed contin- 
uous phase transition is thus that the symmetry group of 
one phase (B in the example) is a sub-group of the other 
(phase A)m& Further restrictions are implied by detailed 
consideration of the LGW expansion (e.g., presence of 
cubic terms, etc.), as is standard^ 

Recent work on related but distinct quantum phase 
transitions has provided explicit theoretical examples of 
non-LGW critical theories^ demonstrating that the vi- 
olation of this conventional "criterion" is more than a 
formal possibility. Unfortunately, there is at present no 
general prescription to supplant the LGW criterion, so 
we are left in the uncomfortable position of being unable 
to solidly argue for or against the possibility of a con- 
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FIG. 6: (Color online) Global phase diagram of the CCM 
model in scenario II of Sec. IV51 The continuous 0(3) tran- 
sition between the Neel* and P-VBS phases is indicated by a 
dashed vertical line. Other notations are as in Fig. |S] 

tinuous quantum critical point (QCP). Instead, we will 
content ourselves here with the LGW analysis. 

It is straightforward to conclude that a continuous 
transition between the CD- VBS and P-VBS states is 
prohibited by the LGW criterion. This can be seen by 
the two lattice reflections 7?-i.2, which map the crossed- 
chains lattice onto itself, i.e., symmetries of the Hamil- 
tonian. Here IZi is the reflection with respect to a hori- 
zontal chain [it corresponds to a link parity operation Pl 
Eq. (|18|l on all vertical chains], and IZ2 is the reflection 
with respect to a horizontal line passing through the cen- 
ters of empty plaquettes [this is a site parity P s Eq. i|17|) 
from the point of view of vertical chains] ; similar reflec- 
tions with respect to vertical lines are accounted for by 
the 7r/2 rotational symmetry (about chain crossings) of 
the lattice. 

Both phases are two-fold degenerate, as can been seen 
from Fig. [S] (and hence can be described by Ising order 
parameters). Their symmetries are distinct. In partic- 
ular, note first that 7£i is a symmetry of the CD- VBS 
phase, but not the P-VBS (it interchanges the two P- 
VBS ground states). Thus the symmetry group of the 
CD- VBS phase is not a subgroup of that of the P-VBS 
phase. Second, 72.2 is a symmetry of the P-VBS phase, 
but not the CD- VBS phase. Thus the symmetry group of 
the P-VBS phase is not a subgroup of the CD- VBS state. 
Since neither symmetry group is a subgroup of the other, 
a continuous LGW transition between the two states is 
not possible, as promised. 

The simplest alternative is a first order transition be- 
tween the two phases, which is always possible, and may 
perhaps be likely. Another possibility, is that, between 
the two states, there is a finite range of coexistence of P- 
VBS and CD- VBS order. Such a coexistence phase can 
have continuous LGW (Ising) transitions to both the P- 
VBS and CD- VBS states. The latter scenario is only one 
of a multitude of conceptually possible phase structures, 
for which we have no physical motivation. We indicate 
this uncertainty by the question mark in Fig. [5] 



B. Scenario II: CD to P-VBS via an intermediate 
ordered phase 

A quite different scenario is suggested by the 
quadrumer-boson approximation of Sec. IIVI an intcrme- 
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diate magnetically ordered phase between the CD-VBS 
and P-VBS states. It was found there that the P-VBS 
state becomes unstable at (J X /J) C 2 = 8/11 as the Jx/J 
ratio is reduced below the planar pyrochlore value of 1. 
The resulting state, depicted in Fig. 0] and Fig. pos- 
sesses long-range magnetic order. We denote it as the 
Neel* state in the following. This magnetically ordered 
state was previously found in the largc-5 approach 3 -, 
where it arises as a result of a quantum "ordcr-from- 
disorder" selection among the large family of degenerate 
(at S = oo) collincar ordered states. The fact that it also 
appears as a result of the triplet softening of the S =1/2 
P-VBS phase of Sec. II VI gives a strong independent argu- 
ment in favor of its stability. Of course, as in the previous 
scenario, the ultimate fate of the Neel* phase is to be de- 
cided by accurate numerical investigations, and here we 
just assume that this ordered phase is indeed the ground 
state in a finite Jx/J window within the (0, 1) interval. 

The transition between the P-VBS and the Neel* 
phases manifests itself as a triplet condensation transi- 
tion within the consideration of Sec. IIVI Since such a 
"soft-mode" transition is the general physical interpre- 
tation of LGW theory, it should not come as a surprise 
that the Neel* to P-VBS transition satisfies the LGW 
criterion for a continuous critical point. In particular, 
the symmetry group of the Neel* state is generated by: 
(1) spin rotations about the ordering axis, (2) transla- 
tions along a diagonal [(1, 1) or (1, —1) in the coordinate 
system of Fig. E] composed with a time-reversal which 
inverts the spins S r — ► —S r , (3) reflections IZ2, and (4) 
7r/2 rotations about the center of a plaqucttc contain- 
ing four parallel spins. Plainly, from Fig. H3 every one of 
these operations is also a symmetry of the P-VBS state; 
hence the symmetry group of the Neel* state is a sub- 
group of that of the P-VBS state. Moreover, the triplet 
condensation amplitude can be identified with the Neel* 
order parameter: an 0(3) vector specifying the direc- 
tion of spin orientation at some reference site in Fig. 
Indeed, an LGW expansion could be developed in this 
order parameter, but we content ourselves with the ex- 
pectation that the P-VBS to Neel* transition is likely in 
the continuous 0(3) universality class. 

Clearly, the Neel* phase cannot survive down to the 
J x /J = point which describes decoupled S = 1/2 
spin chains, as quantum spin fluctuations destroy antifer- 
romagnetic Neel long-range order in individual chains. 3 
The calculations of Sec. II I II demonstrate that at small 
Jx/J, the fluctuating dimerization field e takes over the 
(quasi-classical) spin fluctuations and drives the chains 
into the CD-VBS phase. 

We can again apply the LGW criterion to ask whether 
a continuous CD-VBS to Neel* transition is possible. 
Clearly, the symmetry group of the CD-VBS state can- 
not be a subgroup of the Neel* phase, since the former is 
spin-rotationally invariant. However, the Neel* phase is 
invariant under IZ2 > which, as we saw in the previous sub- 
section, is not a symmetry of the CD-VBS phase. Thus, 
the symmetry group of the Neel* state is not a subgroup 



of that of the CD-VBS phase, and a continuous LGW 
transition between these two phases is prohibited. 

The transition then is likely cither a first-order one 
or proceeds via an intermediate ordered and bond- 
modulated coexistence phase. Such a phase is easy to 
imagine — start with the Neel* state and modulate 
slightly spin exchange couplings along horizontal and 
vertical chains so that "strong" bonds repeat the pat- 
tern of dimers in the CD-VBS phase. This state clearly 
breaks IZ2 , but does preserve the long-range magnetic or- 
der (for sufficiently weak modulation) of the Neel* phase. 
The transition from the CD-VBS to such a "modulated" 
Neel* one is then continuous 0(3) spin-symmetry break- 
ing, while at some higher Jx/J the bond modulation goes 
away via a continuous Z2 transition and one obtains the 
pure Neel* phase. Which of the two possibilities, first- 
order or coexistence, is realized cannot be decided within 
our analytical approach, and for this reason the transi- 
tion between CD-VBS and Neel* phases is marked by a 
question mark in Fig. 



C. Plaquette VBS to Neel transition 

Regardless of the phase structure for J x /J < 1, we 
expect a transition at larger J x from the P-VBS state to 
the conventional Neel state. As is the case with the CD- 
VBS to Neel* transition discussed above, a continuous 
P-VBS to Neel QCP is forbidden within LGW theory (by 
similar arguments, which we refrain from giving explicitly 
for brevity). In this particular case, however, an alert 
reader may note that the two phases appear to be very 
similar to those recently argued to be connected by a 
continuous but non-LGW continuous phase transition, 
deemed a DeconGncd Quantum Critical Point (DQCP), 3 4 

The analogy, however, is not complete. Significantly, 
the checkerboard lattice differs in detail from the square 
lattice discussed in Ref. |^ in its point symmetry group. 
In particular, a tt/2 rotation about a site of the lattice, 
a symmetry of the square lattice, is not a symmetry op- 
eration of the checkerboard lattice. We believe this sym- 
metry distinction is sufficient to destabilize the putative 
DQCP. It is beyond the scope of this paper to fully reca- 
pitulate the arguments of Ref. 0, which would be neces- 
sary to explain this conclusion in a stand-alone fashion. 
Instead, we will sketch these arguments, assuming the 
reader will refer to Ref. for further details. 

The crucial, indeed defining, property of the DQCP is 
an emergent topological conservation law, exactly main- 
tained at the critical fixed point. Specifically, "skyrmion 
number" is conserved by the fixed point theory. This 
is not true microscopically at the lattice level, but is an 
emergent feature of critical theory, as argued in Ref. 
A crucial step in that argument is the remarkable identi- 
fication (due to Haldane2*i) of the skyrmion creation op- 
erator with the columnar/plaquette VBS order param- 
eter. These can be defined through a complex scalar 
field ^ (see Ref. I34F) . We emphasize that although has 



the symmetries of the rather "conventional" VBS order 
parameter, it is not to be viewed as a nearly-free field 
in the LGW sense, but rather some "composite" oper- 
ator in the critical field theory. Under the tt/2 site ro- 
tation above, one finds \& — * and consequently, for 
the square lattice, the skyrmion creation operator can 
appear only in the fourth order in the continuum field 
theory action for the square lattice antiferromagnet, i.e., 
as a perturbation of the form S' — JdrcPx A4RC 'F 4 . A 
variety of arguments indicate then, because of the large 
(fourth) power of \& which appears, A4 is irrelevant at the 
DQCP. On the checkerboard lattice, lacking this tt/2 site 
rotation, the remaining symmetries of the system are in- 
sufficient to rule out the much more relevant "quadratic" 
term S' = JdrcPx A2Re 'F 2 . The presence of the non-zero 
A2 term is, incidentally, also tied to the only two-fold de- 
generacy of the P-VBS state, compared to the four-fold 
degeneracy of the columnar and plaquettc VBS states on 
the square lattice. 

A necessary and sufficient condition for the stability 
of this DQCP on the checkerboard lattice is thus the ir- 
relevance of A2. Ultimately, this can be decided only by 
detailed numerical calculations of the scaling dimension 
of the \I/ 2 operator. Unpublished numerical results^ for 
the easy-plane deformation of the theory (which is more 
numerically tractable) suggest that it is in fact relevant. 
If this conclusion is, as we suspect, true for the SU(2) 
symmetric model, the DQCP is not stable on the checker- 
board lattice. Thus we are led to conclude that there is 
no viable candidate theory for a continuous Neel to P- 
VBS quantum critical point in this model, and that such 
a transition is quite unlikely. 

The simple harmonic analysis of the previous section 
predicts another "soft mode" transition in the 0(3) uni- 
versality class out of the quadrumerized VBS state to one 
with Neel order at J/ J x < 1. The resulting magnetically 
ordered phase is, by its very construction, a coexistence 
region with both P-VBS and Neel order, that is, with 
less symmetry than either phase. This is built into the 
quadrumer boson expansion because all excitations are 
constructed about a background that explicitly has the 
reduced symmetry of the P-VBS state, and there is no 
mechanism to restore the full point group symmetries of 
the checkerboard lattice. Thus we believe the alternative 
possibility of a direct first-order (since the DQCP the- 
ory is unstable) transition from the P-VBS to the true 
Neel state should not be ruled out as a possibility. The 
possible existence of a coexistence region between VBS 
and Neel orders in various models is still a subject of 
some contention. It has been discussed in great detail 
in Ref. l3?l An exact diagonalization study of the quan- 
tum checkerboard antiferromagnet has concluded that, 
if present at all, the co-existence phase is very narrow. 
Clearly more detailed studies of this interesting question 
arc needed. At present, we can only reiterate that a sin- 
gle continuous transition is highly unlikely in view of the 
arguments presented above. 
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FIG. 7: (Color online) Three-dimensional dimerization pat- 
tern of a generalized "quasi-one-dimensional" pyrochlore an- 
tiferromagnet. Also shown, by a long blue arrow, is its "two- 
dimensional" projection, which coincides with the crossed- 
dimer order of Fig. 

VI. BACK TO THREE DIMENSIONS 

Consider now the 3D pyrochlore. Although all bonds 
of a tetrahedron are equivalent by symmetry, it makes 
sense to ask, in analogy with the 2D lattice that we an- 
alyzed in this paper, what would happen if some bonds 
were stronger than others. The particular generaliza- 
tion motivated by the present study of the 2D projected 
model (which can be thought of as a "shadow" of the 
3D pyrochlore on a 2D plane), involves a modified model 
in which two opposite bonds of tetrahedron are strong 
(J) whereas four remaining bonds, connecting the strong 
ones, are weak (Jx). Then, in the limit J x /J "C 1, 
one is back to the problem of strong spin chains cou- 
pled by weak and frustrated inter-chain J x . Now, how- 
ever, chains are arranged in layers: chains are parallel 
to each other (oriented along either x or y direction) in 
each layer, but are orthogonal to those in the layers right 
above and below. That is, chains form a stack of the type 
x — y — x — y... along the vertical (z) direction. Chains 
in one layer do not interact with each other, J x couples 
orthogonal chains from neighboring layers. This is just 
a 3D generalization of 2D situation analyzed in this pa- 
per. It does not introduce any new features, and hence 
the answer for the ground state is straightforward - it is 
spontaneously dimerized into the pattern shown in Fig.0 

Such a generalization is not unrealistic. It appears that 
S =1/2 pyrochlore material GeCu204 has exactly such 
a quasi-one-dimensional structure^ thanks to a strong 
Jahn- Teller elongation of CuOg octahedra along the crys- 
tal c direction^ From the high-temperature tail of the 
uniform spin susceptibility one can estimate the ratio of 
(frustrated) inter-chain exchange J x to the intra-chain 
one J as J x /J ~ 0.16^& At lower temperatures the uni- 
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form spin susceptibility follows that of the spin chain 
down to T c = 33 K, where a small discontinuity is ob- 
served. The specific heat shows a sharp peak at the same 
temperature, suggesting a first-order transition to the 
magnetically ordered state of yet unknown structure (as 
far as we know, Ref . |H is the only experimental study of 
this interesting material at the moment). The theory pre- 
sented in this paper predicts that for a sufficiently small 
J x /J ratio, the ordered state will be replaced by the 
quantum-disordered valence-bond solid shown in Fig. [7] 
Another very interesting realization of "one- 
dimensionality" in the 3D setting is provided by 
the 5=1 pyrochlorc material Znv^O/ijSi There spin 
chains are formed below the structural orbital- ordering 
transition at T c \ = 50 K, as observed in the recent 
neutron scattering experiment^ This is followed by 
the second, magnetic, transition at T C 2 = 40 K. The 
resulting collincar magnetic order can be described as 
follows: 5 = 1 spins order antiferromagnetically along 
the strong (chain) directions in such a way that spins 
along weaker ( J x ) bonds form a "4-spin" pattern^ 
"up-up-down-down-up-up-. . . ." This 3D structure is, 
in fact, rather similar to the 2D one, Fig. 0] found in 
Sec. IV Bl observe that spins on the inter-chain bonds in 
that Figure follow the same "4-spin" pattern of two up 
and two down spins. This is not a coincidence - in both 
cases the classical order can be understood in terms of 
the well-known "order-from-disorder" phenomenon^ 
induced by either quantum^ or thermal fluctuations. 
This analogy suggests that magnetically ordered state 
of GeCu2C>4, observed in Ref. |33, should be similar to 
that in the low-temperature phase of ZnV2C>4 — clearly 
more experimental and theoretical studies of this ques- 
tion are required. The analogy does not apply in the 
Jx/ J - * limit where decoupled 5 = 1 chains, although 
in the quantum-disordered phase with a finite spin gap 
A ~ 0.4 do not break translational symmetry. This 
is in contrast with 5 = 1/2 chains studied in this paper; 
the decoupled limit is characterized by the crossed-dimcr 
order, Fig. [7| which does break translational symme- 
try. Properties of the 3D phase transitions between 
quantum-disordered and ordered phases constitute 
another interesting theoretical problem which we leave 
for future studies. 



VII. CONCLUSIONS 

The main result of this work consists in the predic- 
tion of the novel VBS phase, the crossed-dimer phase, 
illustrated in Fig. [3] and Fig. □ This new VBS phase 
arises in the ID limit of the model as a result of the 
frustration-fostered competition between classical (repre- 
sented by the staggered magnetization N) and quantum 
(represented by the staggered dimerization e) ordering 
tendencies. Our analysis is based on the careful pcrturba- 
tive implementation of the well-known SU(2) symmetry 
of the g matrix field of the WZW model, which provides 



rigorous field-theoretical description of the low-energy 
sector of the 5 = 1/2 isotropic Heisenberg chain. This 
symmetry is made explicit by the OPE Eqs. l|2"5|l - l|2"8j l. 
which show transformation properties of the low-energy 
"quantum triad" { Jr/l, N, e}. As shown in Sec. lIIII con- 
sistent implementation of these OPEs requires a careful 
treatment of the often- neglected gradient terms ("non- 
primary fields" in the conformal field theory nomencla- 
ture, such as d x N) which link together quantum fluctua- 
tions of the conserved spin current with those of the stag- 
gered magnetization and dimerization fields. Once this 
is done, straightforward inter-chain perturbation theory 
leads to the frustration-generated interaction of dimer- 
izations on the crossing chains, Eq. (|40|l . 

Our finding of the CD- VBS phase in the J x < J limit 
of the checkerboard antifcrromagnct eliminates a previ- 
ously proposed^ sliding Luttingcr liquid phase as the can- 
didate for the ground state. Like many others, that work 8 
overlooked the crucial role of the gradient terms in the 
analysis of the frustrated inter-chain interaction between 
critical 5 = 1/2 Heisenberg chains. 

It is also worth pointing out here that our calculation 
clarifies previous somewhat inconclusive "sightings" of 
the decoupled-chains phased that arise in a widely-used 
large-A^ approach^ to frustrated spin models. By its 
very construction, that technique fails to account, at the 
leading N = oo level, for the fluctuation-generated resid- 
ual dimer-dimer interaction in the anisotropic ID limit 
(although one does expect finite 1/N corrections to the 
inter-dimer interaction to appear once the fluctuations of 
the compact gauge field are accounted for A) 

We have also presented a global phase diagram of the 
CCM (Sec. |Vj. Although phcnomcnological in nature, 
our analysis stresses the importance of lattice symmetries 
in delineating the order of possible direct transitions be- 
tween various quantum (CD- and P-VBS) and classical 
(ordered Neel and Neel* ) phases of the CCM found in this 
and previous studies. We find that most of such transi- 
tions are required to be of the first-order type, or proceed 
via an intermediate co-existence phases, as illustrated in 
Figs.|Sland|SJ This claim concerns even the relatively well 
studied P-VBS to Neel transitioniiiiSii^ and clearly calls 
for more numerical (as well as analytical) investigations 
of this interesting question. 

Last but not least, we have also presented a sim- 
ple but intriguing extension of the approach to the 
anisotropic three-dimensional pyrochlore antiferromag- 
net, which may be relevant to both 5=1/2 and 5=1 
pyrochlore-based magnetic materials iSiial We hope that 
this interesting connection will inspire new experiments 
in this exciting area. 
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APPENDIX A: CORRELATION FUNCTIONS 
1. Fermionic formulation 

Formally, the WZW model in Eqs. © an d (O de- 
scribes all the low-energy properties of the S = 1/2 
Heisenberg spin chains, and is the starting point for a 
perturbative analysis of interchain coupling terms. Un- 
fortunately, as remarked in Sec. Ill Bl this action and 
Hamiltonian, while compact and self-contained, are not 
directly useful for concrete calculations. Indeed, in two- 
dimensional conformal field theory, analysis often pro- 
ceeds without any action or Hamiltonian, solely on the 
basis of algebraic operator product relations that essen- 
tially specify the field content and correlation functions of 
the theory. We will ultimately proceed by using operator 
product relations along these lines. In the crossed chains 
model, however, due to the local nature of interactions 
between perpendicular chains, we need to pay particular 
attention to proper short distance regularization of the 
theory. It is therefore desirable to have a more concrete 
formulation of the theory. 

To do so, we use the well-known phenomena of spin- 
charge separation in one-dimensional spin- 1/2 Dirac 
fcrmions. In particular, using bosonization, one may 
show that all degrees of freedom (all states and opera- 
tors) of such Dirac fermions are identical to those of a 
free scalar "charge" boson and a spin sector described by 
the WZW SU (2)i theory of interest to us. Moreover, the 
free Dirac fermion Hamiltonian, 

Hd = iv J dx (*l t ,d x * Li , - <A*fl,.) . ( A1 ) 

can be expressed as the sum of decoupled spin and charge 
Hamiltonians, Hd = Hwzw + H p , with 

H P = lJ dx [( d *<P P ) 2 + (W 2 ] . ( A2 ) 

where (p p ,i9 p are conjugate "charge" boson fields 
[ip p (x),'dp(x')] = —isign(x—x')/2. The right-moving and 



left-moving fermion operators may, if desired, themselves 
be re-expressed in terms of the gapless spin and charge 
fluctuations, which is commonly and conveniently done 
with the help of abelian bosonizationiS^ We do not, how- 
ever, require these expressions here. 

Because the spin degrees of freedom in Hwzw (de- 
scribed by the S[/(2)-valucd field g) are independent of 
the charge boson fields, we may safely replace Hwzw by 
Hd at the price of enlarging the physical Hilbert space 
to include these charge degrees of freedom, without how- 
ever affecting the spin physics in any way. In particular, 
regardless of the spin interactions we add to Hwzw, the 
charge sector remains always in its ground state, the vac- 
uum of Eq. (|A2j) . 

The fermionic formulation provides a convenient way 
to calculate while simultaneously regularizing the short- 
distance properties of the theory. This is accomplished by 
representing the important operators in the spin sector in 
terms of the Dirac fermions. First, the spin currents may 
be fully and simply re-expressed in terms of the Dirac 
fermions: 

4 = *L%^> ■>> "iM*! - (A3) 

We can thereby construct the uniform magnetization, 
J = Jr + Jl in terms of fermions. Observe that J is 
invariant under chiral U(l) charge symmetry 

*B/L,.-e < **/- 6 * Ji/£ , a (A4) 

with independent phases \R XL- This invariance im- 
plies independence of J from the charge sector, Eq. I|A2(1 . 
no matter what its precise form is. 

The second important operator in the spin sector is 
the staggered magnetization. Unlike the uniform mag- 
netization, however, it does not have a simple expres- 
sion in terms of the fermions. Instead, one may define 
a fermionic "staggered" (corresponding to the 2kp com- 
ponent for physical electrons in the Tomonaga-Luttinger 
model) spin density: 

^■ = *k.% 1 *^ + *L.% 1 *^- (AS) 

Using standard bosonization relations, one finds 

N F = N cos V2ttlp p . (A6) 

Hence the fermionic staggered spin density, Eq. I|A5|I . 
reproduces the desired staggered magnetization, iV, but 
multiplied by a factor involving the charge boson (p„. 

Observe that Np is not invariant under the chiral 
charge symmetry Eq. (|A4|) and this implies that it does 
couple to the charge sector of Eq. (|A1|) . as is explicitly 
stated in Eq. (|A6|) above. (It is invariant under the di- 
agonal subgroup of Eq. (|A4|) with %R = XL which only 
reflects that it conserves the total charge of the system, 
but not the charges of right- and left- moving chiral sec- 
tors independently) 
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Similarly, wc may define a fermionic staggered dimer- 
ization operator, 



by 



;(*L 



Bosonization gives 



ep = e cos \/2Tzip p 



(A7) 



(A8) 



Like for the staggered magnetization, ep defined from the 
Dirac fermions reproduces the staggered dimerization op- 
erator e in the spin sector, but with an undesired multi- 
plicative charge factor. Here we observe again that this 
multiplicative charge factor follows from the fact that ep 
does not respect Eq. (|A4|) . similarly to Np above. 

We will see below that, for our purposes this second 
factor in Eqs. i|A6(l and i|A8() is innocuous, so that we 
may use Np,ep in place of the true N,e fields. We 
stress this is not generally true, and such a replacement 
is possible only under special circumstances we outline 
below. For example, because Np, ep contain charge fluc- 
tuations, their correlation functions are not the same as 
those of the N, e operators, and hence the latter can- 
not be evaluated using free fcrmion apparatus. However, 
the fermionic operators have the same operator products 
with the spin current field J as the corresponding spin 
operators. This is because the spin current, as discussed 
above, does not depend on charge sector, and as a result 
the "charge pieces" of Np and ep stay "inert" during 
OPE calculations as outlined in Appendix [51 As it turns 
out, such operator products are all we need in order to 
solve the problem. Thus, we can evaluate these crucial 
operator products using those of the free fermions. Fur- 
thermore, we can regularize the problem by adopting a 
short-distance cutoff for the Dirac fermions. This is par- 
ticularly convenient, as a simple ultraviolet regularization 
is available in the Dirac formulation which keeps SU(2) 
spin rotation symmetry manifest at all stages. 



2. Fermion Green's function. 

To proceed with our program we need then the Green's 
function of free fermions ^ r /l ■ As the inter-chain inter- 
action, Eq. \23$ . is the sum of local terms over crossings 
of the lattice, one needs to be careful about the high- 
energy regularization of the continuum theory. We have 
found that the most natural and physically appealing reg- 
ularization is provided by the "soft-momentum cutoff" 
scheme4£*2ii In that scheme cutoff a is introduced via 
e -a\q\ f ac t or i n the mode expansion of collective bosonic 
fields. In going from the momentum to the coordinate 
space, cutoff dependence on a transforms into that on 
q = a/v in the temporal (r) direction. Adopting this 
scheme, one finds that the zero temperature (T = 0) 
Green's function for the right-moving fermions on the 
same chain (chain index is omitted for brevity) is given 



F r (x,t) = -(f* RtS (x,T)*+ tS (0,0)) 
1 



2ttv(t — ix/v + aa T ) 



(A9) 



where a = a/v is the ultraviolet cutoff in temporal di- 
rection and ay = sign(r) = O(r) — 0(— r) is the sign- 
function. For the left-moving ones 



F l (x,t) = F r (x,t) = F R (-x,r) = -F r (x,-t) 



1 



(A10) 



2ttv(t + ix/v + acr T ) 
It is often convenient to introduce complex coordinates 
z = t + ix/v, z = r~ix/v (All) 
in terms of which 



Fr{z) 



1 



2irvz 



, F L {z) 



2irvz 



(A12) 



It is important to keep in mind, however, that expressions 
of Eqs. (|A9|) and (|A10|) are "superior" to Eq. I|A12|) in 
that they contain an explicit cutoff prescription. 



APPENDIX B: DERIVATION OF OPE EQ. (|2?j) 

We need Eqs. 1|A3|) and l)A5jl and the single-particle 
Green's function, Eqs. I|A9|I and (|A10|) . Consider the 
product of J^,(.ti,ti) and N F (x2,T2) at nearby points 
(such that \zi — z% | ~ a with = Tj + ixi/v). To shorten 
notations, wc denote (xi,Ti) = 1, etc. Now use Wick's 
theorem for free fermions, as appropriate for Eq. i|Alfl . 
as well as the identity a a a h = 5 ab + ie abc a c , to reduce the 
product of four ^'s to that of two of them and Fr 

r R {l)N F {2) = -\F R {l-2) 

x^e abc alA^J^L.A2) + ^t s (l)^ R , s/ m 
+ 5 o6 [*+ s (l)* i)S (2) - ^(l)**,.^)]) , (Bl) 

where summation over repeated spin indices s and s' is 
assumed. Fusing (xi,Ti) and (22, T2) points, one finds 
that the coefficient of e abc is just 2N F (2), whereas that 
of S ab is — 2iep(2), the staggered dimerization operator; 
see Eq. (|A7I) . Hence 



J a R {l)N b F (2) 



ie abc N c F {2) - iS ab e F (2) 
4ttv(zi - z 2 + ao Tl - T2 ) 



(B2) 



The second OPE, between Jl and Np, is obtained by 
replacing Fr — > Fl and changing the sign of the last 
term in Eq. I|B2|) as is readily verified by the explicit 
calculation. Thus 



Jl(l)N F (2) 



ie abc N c F {2) + iS ab e F {2) 
4ttv(zi - z 2 + aa Tl - T2 ) 



(B3) 
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Observe now that the multiplicative charge factor 
cos y/2TTip p from Eqs. (|A6() and (|A8() appears on both 
sides of Eqs. I|B2|) and l|B3|) above. Dividing by it, we 
obtain the OPE quoted in Eq. |J27|) 



J%(z)N b (w,w) = 
J a L {z)N\w,w) = 



ie abc N c {w,w) - iS ab e{w,w) 

At:v(z — w) 
ie abc N c (w, w) + iS ab e(w, w) 

47T«(z — w) 



(B4) 



It is instructive to repeat the same calculation us- 
ing abclian bosonization. Consider, for example, 
J%(x\, Ti)N x (x2, T2). Abclian bosonization tells us that 

1 



J£(l) = -=0 xi * fl)ir (l), (B5) 

JV«(2) = _^_f e -*^I*«.»(i)-»i.-( 1 )]+H.c.V (B6) 
47ra V / 



where <&r/l,<t are the spin components of chiral right/left 
bosons, and the scale factor A = (cos v / 27r<Pp)p is the ex- 
pectation value of the charge field. The brackets (...) p 
denote average with respect to the charge Hamiltonian 
H' p which includes now four-fermion umklapp term re- 
sponsible for the opening of the charge gap (as a result 
of which A ^ fl), see Appendix A of Ref. [23 for more 
details. 



Thus charge fluctuations are explicitly separated from 
spin ones in the manipulations to follow. To fuse Eq. 1B5I) 
and Eq. (|B6(I . we need 



i / 3* H , <T (2) 



(B7) 

which is obtained by expanding exponential on the left- 
hand side of Eq. (|B 7|) . pairing <6r i<t (1) with <&r i<t (2) in 
all possible ways (Wick's theorem for free chiral spin 
bosons), and collecting the rest of the series back into 
exponential. Using Eq. (|B7|I and the fact that right and 
left bosons are independent from each other, we find 

J Z R (1)N X (2) = d Xl (f$RAl)®R.A2))N y (2). (B8) 

Finally, 

(f$R )CT (l)$ fli(7 (2)) = --/-ln[a + a T1 _ T2 (l 1 -f 2 )] (B9) 

47T 



for free chiral bosons. Hence 



r R (i)N*(2) 



iN y {2) 



4irv(zi — z 2 + a<r Tl - T2 ) ' 



(BIO) 



which is just one component of Eq. (|B2|) . To get 5 ab term 
of Eq. (|B2|I one has to consider explicitly OPE of, say, 
J z and N z fields. 
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